 #include "hamilton.hpp"

hamilton::hamilton(){}

hamilton::hamilton(matrix<double> H){
  H_=H;
}

hamilton::~hamilton(){}

matrix<double>  hamilton::operator()(){
  return H_;
}

matrix<double>  hamilton::get_op(){
  return H_;
}

hamilton::hamilton(const hamilton& h){
  if(this != &h){
    H_=h.H_;
  }
}
hamilton& hamilton::operator=(const hamilton& h){
  if(this != &h){
    H_=h.H_;
  }
  return *this;
}
// vector<std::complex<double> > hamilton::operator*(matrix<std::complex<double> > h,
// 					vector<std::complex<double> > psi){

// }
